Double- Well Quantum Tunneling Visualized via Wigner's Function 
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We investigate quantum tunneling in smooth symmetric and asymmetric double-well potentials. 
Exact solutions for the ground and first excited states are used to study the dynamics. We introduce 
Wigner's quasi-probability distribution function to highlight and visualize the non-classical nature 
of spatial correlations arising in tunneling. 
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I. INTRODUCTION 

Quantum tunneling was first discussed by Fricdrich 
Hund in 1927 when he considered the ground state of 
a double- well potential 0, Quantum tunneling is a 
microscopic phenomenon where a particle can penetrate 
into and in some cases pass through a potential barrier, 
although the barrier is energetically higher than the ki- 
netic energy of the particle. This motion amounts to 
particles penetrating areas where they have "negative ki- 
netic energy" and is not allowed by the laws of classical 
dynamics Q. Double- well potentials can be used for the 
study of tunneling, as the central hump separating the 
two wells constitutes a tunneling barrier, for eigenstates 
lower in energy than the maximum height of the barrier. 
They are also used for the study of molecular structures, 
for example in the ammonia molecule Q . The concept of 
quantum tunneling is central to the operation of scanning 
tunneling microscopes [H, [f| . 

In this paper we investigate tunneling in partially ex- 
actly solvable double-well potentials Q through the use 
of the Wigncr quasi-probability distribution function. 
We model tunneling based on the dynamics of the low- 
est two wave functions, the quantum mechanical ground 
state and the first excited state 0- We also present the 
probability distributions with respect to momentum and 
position, and illustrate the presence of quantum coher- 
ences, which give rise to interference fringes with negative 
values of Wigner's quasi-probability distribution. 

In the next section, we provide a self-consistent dis- 
cussion of the partially exactly solvable symmetric and 
asymmetric double- well potentials, first established in 
In section IIIII we introduce Wigner's quasi-probability 
distribution function and illustrate, for each case of the 
double-well potential, time-evolution of the associated 
Wigner quasi-probability distribution, and probability 
distributions with respect to position and momentum. 
Subsection MI Al considers negative regions of Wigner's 
quasi-probability distribution function and allows us to 
visualize non-classical spatial coherences in tunneling, fi- 
nally we conclude. 



II. PARTIALLY EXACTLY SOLVABLE 
DOUBLE- WELL POTENTIALS 

A. The Schrodinger Equation 

A family of partially exactly solvable double-well po- 
tentials was introduced in 1995 by Caticha Q; for these 
the ground and first excited states (and in some cases 
states energetically above those) can be computed ana- 
lytically in simple closed forms. Note that despite consid- 
erable efforts.only few fully solvable smooth potentials 
arc known [8, 9]; no fully solvable smooth double- well po- 
tential has yet been found. Such partially solvable mod- 
els may well prove to be very useful as benchmarks for 
computer code and numerical tests. 

The Schrodinger equation (for the two lowest states, 
n = and n = 1) is 

h 2 

—i; n "{x) + {E n -V(x)]ip n (x) = 0. (1) 
Zm 

Throughout this paper we choose H 2 /2m = 1. Similarly 
to super-symmetric quantum mechanics [1, Q, where a 
superpotential is defined in order to compute the ground 
state and the corresponding potential, in Q a multiplier- 
function (f> is defined, which relates the ground and first 
excited states by 

V-i = Ho • (2) 
Substituting Eq. into Eq. JT]) for n = 1 yields 

(#o)" + [Ei -V]Hi Q = 0. (3) 
Eq. (JTJ) for n = multiplied by (f> is 

Ho + [E - V] M'o = . (4) 

Upon subtraction of Eq. (UJ from Eq. © , one obtains a 
function \ 



AE<f) 



2c5' 



din ipo 
dx 



(5) 



where AE = E\ — Eq is the energy splitting between 
ground and first excited states. Eq. ([S]) resembles a su- 
perpotential in the context of super-symmetric quantum 
systems [1, Q . The corresponding ground state tpo there- 
fore is 



* D.Kakofcngitis@hcrts.ac.uk 
t O.Stcucrnagcl@hcrts.ac.uk 



ip (x) = AAexp J x(x')dx'^ , 



(0) 



2 



here N is a normalization constant. Rearranging Eq. §4§ 
determines the double- well potential V, up to an additive 
constant 



En — 



o = X - X 



This is illustrated in Fig. [TJ 



En 



(7) 



and the first excited state ipi is 



V>i(0)(e 



l)(a-b) 



hi- 



( e 2(b+a)x + l)(a-b) + (a + b){e 2ax + e 2bx ) 



(13) 

The two wave functions have odd and even parity, see 
Fig. [2 





FIG. 1. The double-well potential V (solid line), the func- 
tion x (dashed line) and the multiplier-function <j> (dotted 
line). For the symmetric case, (a, cf> = sinh(aa-)/ cosh(fo:r), for 
the asymmetric case, (b, (f> — a + tanh(/3a;). The symmetric 
case with parameters a = \J — En, b = \/—Ei, Eq = — f and 
Ei = —0.999, features odd x an( f <t> functions, whereas for 
the asymmetric case (here with a = 0.9, /3 — 1, En = and 
Ei = 1) these symmetries are lost. 



B. Symmetric Double- Well Potential 

For the case of a symmetric double-well potential, as 
displayed in Figs. [T] a and [H we use the multiplier- 
function 



<t>(x) 

where (with E < E 1 < 0) 



sinh(ax) 
cosh(6a;) 



a = \/—Eq and b = sf—E\ 
The function x then takes the form 



(8) 



(9) 



sinh(ax) (2b 2 - 2a 2 cosh 2 (6x)) 



6sinh(aa;) sinh(26a;) — 2acosh(aa;) cosh (bx) 
abcosYi(ax) sinh(26;r) 



6sinh(aa;) sinh(2&ir) — 2acosh(ax) cosh 2 (&x) 
and the corresponding symmetric potential V is 

2 (b 2 - a 2 ) (a 2 cosh 2 (fc.T) + b 2 sinh 2 (aa;)) 



(10) 



V(x) 



(acosh(acc) cosh(6x') — bsmh(ax) sinh(fex)) 2 



The ground state ipo (note typographical error in corre- 
sponding Eq. (19) of ref. Q), is 



tp (x) 



M0)( 



e 2bx + l)(a - b)e a 



( e 2(b+a)x + X )( a -b) + {a + b){e 2ax + e 2bx ) 

(12) 



FIG. 2. (Color Online) The symmetric double-well potential 
V (solid blue line), the corresponding ground state i/jn (dashed 
red line) and the first excited state 4>i (dashed green line), for 
E = -1 and Ei = -0.9. 



C. Asymmetric Double- Well Potential 

With the multiplier-function (f> 

4>{x) — a + tanh(/?a;) 
the function x takes the form 



(14) 



AE 



x (x) = ^ [sinh(2/3x) + 2a cosh 2 (fix)] - (3 tanh(/3x) 

(15) 

and the corresponding asymmetric potential V is 

V(x) = f3 2 - AEa sinh(2/3x) 



cosh 2 (/3x) I — -5-asinh(2/3x) 
V 4 /3 



AE 2 



2AE 



AE 2 3 
+ 4^2" ( a * + l ) cosh 4 (^) + -AS + E {) . (16) 

The corresponding lowest two wave functions are 

ipo(x) = V'o(O) cosh(/3x) 



x exp 



1^ ( cosll2 (' 5:z; ) + a ^ x + \ sinh(2^x) 



(17) 



(11) and 



i\)\{x) = V'i(0)(acosh(/3a;) + sinh(/3x)) 
AE ( , , , „ , „ a 



x exp 



v cosh 2 (/3x) + a/3x + ^ sinh(2/3x^ 



(18) 
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Fig. [3] illustrates that the potential's asymmetry strongly 
modifies the shape of the lowest two wave functions as 
compared to the symmetric case, with a = 0, in Fig. [2] 



and momentum 
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FIG. 3. (Color Online) The asymmetric double-well potential 
V (solid blue line), the corresponding ground state tpo (dashed 
red line) and the first excited state ipi (dashed green line) , for 
a = 0.9, = 1, E = and AE = 1. 

To investigate the tunneling dynamics we use the nor- 
malized superposition state 



ty(x,t) = sin(0)exp 



-iE t 



tpo{x) 



cos(6>) cxp 



-iErt 



M*)> ( 19 ) 



with the weighting angle 9 £ [0, . . . , tt/2]. The energy 
splitting AE, gives rise to the beat period (or recipro- 
cal barrier-tunneling rate HH) T = 2irh/AE. In other 
words, T is the time needed for a quantum particle ini- 
tially localized in, say, the left well, to perform a full 
oscillation: left-right- left. The larger the energy splitting 
the shorter the beat period. 



III. WIGNER'S FUNCTION 

Eugene Wigner introduced his quasi-probability distri- 
bution function W(x,p; t) in 1932 for the study of quan- 
tum corrections to classical statistical mechanics PJ|. It 
is a generating function for all spatial auto-correlation 
functions of a given quantum mechanical wave function 
ip and defined as 0, 0, [TJ 



W{x,p;t) 



1 



+ y, i)*(x - y, t)e * dy 



(20) 

where x and y are position variables and p the momen- 
tum. The Wigner quasi-probability distribution function 
is a real-valued phase-space distribution function, it can 
assume negative values which is why it is referred to as 
a gizasz-probabihty distribution. 

Its marginals are the quantum-mechanical probability 
distributions of position 

/>oo 

P(x,t) = \*(x,t)\* = W(x,p;t)dp (21) 



p(p,t) = i*(p,t)r 



W(x,p;t) dx . (22) 



It is normalized J_ W(x,p;t) dp dx = 1 and 

the overlap of the Wigner functions W$(x,p;t) and 
W x (x,p;t) of two distinct quantum states, ip(x,t) and 
x{x, t), yields the magnitude of their wave function over- 
lap squared [ll[ 



W^(x,p;t)W x (x,p;t) dp dx = — \(ip\x)\ ■ 

(23) 

Fig. 2] shows plots of the time evolution of the Wigner 
functions for symmetric and asymmetric double- well po- 
tentials and the associated marginals P(x, t) and P{p, t); 
all Wigner functions and marginals P(p, t) had to be de- 
termined through numerical integrations. 



A. Negative values of Wigner function 

As can be seen from the projections of the Wigner 
quasi-probability distributions in Fig. [4j a particle that 
exists in two places at the same time shows interference 
fringes in its momentum probability distribution P(p, t). 
These are incompatible with a single humped position 
distribution P(x, t) in conjunction with a positive semi- 
definite phase space probability distribution. A phase 
space distribution simultaneously yielding such marginals 
has to contain regions with negative values. These nega- 
tive regions are indicators of the non-classical character 
of the spatial coherences of the wave functions [l2[ and 
are frequently studied in experiments (l3l |. They arise, 
e.g., whenever the wave function spreads out over both 
wells of the double-well potential. 

To avoid confusion, we would like to emphasize that 
negative regions of the Wigner function, as a generat- 
ing function for spatial auto-correlation functions, rep- 
resent non-classical behavior of spatial correlations but 
not of the non-classical behavior of tunneling associated 
with "negative kinetic energy" . The interference fringes 
of the Wigner function, appearing roughly in the area 
where tunneling occurs, represent the non-classical spa- 
tial coherence of the wave functions located in both wells 
simultaneously. 

Fig. [S] illustrates the changes Wigner's functions un- 
dergo when changing the potential wells from separate 
wells to merged single troughs (displayed in Fig. [B]): 
Since the coupling between the wells increases, so does 
AE. Concurrently the peaks of the spatial wave func- 
tions move together which increases the fringe spacing 
in the associated momentum representation, visible as a 
reduction of the spatial frequency of Wigner's functions' 
interference patterns. 
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FIG. 4. (Color Online) Time evolution of the Wigner quasi- 
probability distribution, for the symmetric double- well poten- 
tial V of Eq. jlT} for E = -1 and Ei = -0.999, shown in 
subfigures a, b and c, and for the asymmetric double-well 
potential V of Eq. |16} for a = 0.9, p = 1, E = and 
-Bi = 0.001, shown in subfigures d, e and f, with weighting 
angle 9 — n/4. The projections of the Wigner function are 
the probability distributions with respect to position P(x,t) 
(black) and momentum P(p,t)/3 (brown - shrunk by "1/3" 
for clarity and visibility). 



IV. CONCLUSION 



We have considered the effects of tunneling in smooth 
partially exactly solvable double- well potentials 0. We 
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FIG. 5. (Color Online - Color scheme identical to Fig.UJs) The 
Wigner quasi-probability distribution for a range of values of 
the energy splitting AE, for the symmetric double- well poten- 
tial V of Eq. jni) for E = -1 and AE = 0.25, 0.5 and 0.75, 
shown in subfigures a, b and c, respectively, and for the asym- 
metric double-well potential V of Eq. ()16|1 for a = 0.9, /3 = 1, 
Eo = and for A_E = 0.5, 4 and 8, shown in subfigures d, e 
and f, respectively, with weighting angle 6 — 7r/4 and t — T/4. 



illustrated the behavior of the associated wave func- 
tions and Wigner's quasi-probability phase-space distri- 
butions. Wigner's functions can assume negative values 
representing non-classical spatial coherences of the wave 
functions, these were shown to arise in the case of tun- 
neling and analyzed and interpreted. 
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FIG. 6. (Color Online) The symmetric double-well potential 
V of Eq. HIl), a, for E = -1 and AE = 0.25,0.5 and 0.75, 
color blue, red and green, respectively, and the asymmetric 
double-well potential V of Eq. (flB). b, for a = 0.9, f3 = 1, 
Eq — and AE — 0.5, 4 and 8, color blue, red and green, 
respectively, and their corresponding ground and first excited 
states, ipo and tpi. 
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